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Q ■ Abstract 

O ' 

O ■ The lifetimes of B + and B mesons are measured using a sample of 

150,000 hadronic Z° decays collected by the SLD experiment at the SLC 
between 1993 and 1995. The analysis identifies the semileptonic decays of 
B mesons and reconstructs the B meson decay length and charge by vertex- 
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ing the lepton with a partially reconstructed D meson. This new method 
results in a sample of 634 (584) charged (neutral) decays with high charge 
purity. The ratio of B + : B° decays in the charged (neutral) sample is 3:1 
(1:3). A maximum likelihood fit yields t b + = 1.6llo;i2(stat) ± 0.07(syst) ps, 
t b0 = 1.56±B;i| ± 0.10 ps, and t b +/t b o = 1.03lg;l| ± 0.09. 
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According to the spectator model of heavy hadron weak decay, the heavy 
quark decays independently of the other quarks in the hadron. Therefore, 
this model predicts that the lifetimes of all hadrons containing a given heavy 
quark Q are equal. However, the hierarchy observed in the charm system, 
r D + > r D + ~ r D o > t a +, indicates the need for corrections to this model. 
Such lifetime differences are predicted to scale with I/tuq |l[ and therefore 
are expected to be less than 10% in the 6-quark system. Measurements of 
the B + and B° lifetimes provide tests of this prediction. Finally, precise 
measurements of exclusive B meson lifetimes are necessary to extract the 
CKM matrix element V c b- 

The measurements presented here use a sample of 150,000 hadronic Z° 
decays collected between 1993 and 1995 by the SLD experiment at the SLC. 
The analysis uses a new technique to identify the B hadron charge by re- 
constructing the charged track topology of both B and cascade D vertices 
in semileptonic B decays. Most previous measurements of the B + and B° 
lifetimes || are based on samples of semileptonic decays in which the lepton 
is identified and the D^ meson is fully reconstructed. The lifetime measure- 
ments then rely on assumptions concerning the B + and B° content of the 
D°Xl + u and D^*'~Xl + u samples. In contrast, the more inclusive technique 
used here only relies on the simple difference of total charge between B + and 
B° decays. However, this technique requires very good vertexing to assign 
tracks correctly to secondary vertices. 

The analysis uses the calorimetry and tracking systems (for details see 
Ref. ||). The Liquid Argon Calorimeter (LAC) is used to reconstruct jets 
from energy clusters and perform electron identification with maximal effi- 
ciency for | cos#| < 0.72. The Warm Iron Calorimeter (WIC) provides effi- 
cient muon identification for | cos#| < 0.60. Tracking is performed with the 
Central Drift Chamber (CDC) and the CCD pixel Vertex Detector (VXD) 
with 94% total reconstruction efficiency for |cos#| < 0.74 (including VXD 
hit linking). The impact parameter resolution for high- momentum tracks is 
measured using Z° — > fi + n~ decays to be 11 /xm in the plane perpendicular 
to the beam axis (r0 plane) and 38 /im in the plane containing the beam 
axis (rz plane). 

The decay length is measured relative to the position of the micron-size 
SLC Interaction Point (IP) which is reconstructed in the rip plane with a 
precision of a r ^ = (7± 2) /im using tracks in sets of ~ 30 sequential hadronic 
Z° decays. The z position of the IP is determined on an event-by-event basis 



with o z ~ 52 /iin for b b events using the median z position of tracks at 
their point-of-closest-approach to the IP in the r<p plane. 

The measurements rely on a Monte Carlo simulation based on the JET- 
SET 7.4 event generator j| and the GEANT 3.21 detector simulation pack- 
age H to determine the charge separation purity and to extract the lifetimes 
from the decay length distributions. The 6-quark fragmentation follows the 
Peterson et al. parameterization ||. B mesons (baryons) are generated with 
mean lifetime r = 1.55 ps (r = 1.10 ps). B meson decays are modelled 
according to the CLEO B decay model 0] tuned to reproduce the spectra 
and multiplicities of leptons, charmed hadrons, pions, kaons, and protons, 
measured at the T(4S) |§. Semileptonic decays follow the ISGW model |§ 
including 23% D** production. B baryon and charmed hadron decays are 
modelled using JETSET with, in the latter case, branching fractions tuned 
to existing measurements [TO . 



The initial step in the event selection is to select electron and muon candi- 
dates using the measured track parameters as well as measurements from the 



LAC and WIC respectively (see Ref. [II] for further details). To enhance the 
fraction of Z° — > b b events with little loss in efficiency, lepton candidates are 
required to have total momentum p > 2 GeV/c and momentum transverse 
to the nearest jet > 0.4 GeV/c (jets are found using the JADE algorithm |12 



with y cut = 0.005). These cuts yield a sample of ~ 34, 000 event hemispheres, 
with an efficiency of ~ 75% for semileptonic B decays within | cos^l < 0.6 
determined from our Monte Carlo simulation. 

The secondary vertex reconstruction proceeds separately for each event 
hemisphere containing a lepton, and uses a multi-pass algorithm that op- 
erates on those tracks that have at least one VXD hit and are not from 
identified 7 conversions, or K® or A decays. Tracks are initially classified 
as primary unless their 3-D impact parameter significance with respect to 
the IP is > 3.5 a and p > 0.8 GeV/c, in which case they are classified as 
secondary. 

In the first pass, the hemisphere containing the lepton candidate is re- 
quired to include no more than four secondary tracks (excluding the lepton) 
and a candidate D vertex is constructed using all such tracks (vertex cuts 
are defined below). The D trajectory, found from the D vertex and the total 
momentum vector of tracks included in the vertex, must intersect the lepton 
to form a valid one-prong B vertex solution. If this is successful, an attempt 
is made to form a two-prong B vertex by attaching one primary track to the 
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lepton near the point of intersection. This first pass identifies 91% of the fi- 
nal candidates. These candidates are allowed to be modified by searching for 
one or two primary tracks that can be added to the existing D vertex. This 
search is successful for 40% of the candidates. In case of multiple solutions, 
we select the one with the largest number of tracks and if more than one still 
remains, we select that with the smallest impact parameter between the D 
trajectory and the lepton or two-prong B vertex. A second pass is performed 
if no first pass candidate is identified. Here, a search is made for solutions in 
which one secondary track makes a valid two-prong B vertex with the lep- 
ton, the remaining secondary tracks form a D vertex, and the D trajectory 
intersects the B vertex. Multiple solutions are handled as described above. 

The requirements to form a D vertex are: the number of tracks is < 4; 
the absolute value of the charge < 1; the mass (charged tracks assumed to be 
it's) < 1.98 GeV/c 2 ; the vertex displacement from the IP > 4cr and < 2.5 cm; 
and the vertex \ 2 (2,3,4 prongs) < (4, 12, 20). The requirements to form a B 
vertex are: the absolute value of the total charge (B+D tracks) < 1; the mass 
> 1.4 GeV/c 2 ; the observed decay length > 0.08 cm and < 2.4 cm; and the 
momentum of the non-lepton track (if any) > 0.4 GeV/c. The requirements 
for the D vertex to be linked to the B vertex are: the signed distance between 
D and B vertices > 200 /im; for one-prong B vertices, the distance of closest 
approach of the D trajectory to the lepton < (130, 100, 70) /im for (2, 3, 
4) prong D vertices; for two-prong B vertices, the three-dimensional impact 
parameter of the D trajectory with respect to the B vertex < 200 /im. 

The algorithm yields 783 charged and 584 neutral semileptonic B decay 
candidates. The topological breakdown is given in Table [l]. The efficiency 
for reconstructing a semileptonic B decay is estimated from the simulation 
to be 24% for decays with an identified lepton within | cos#| < 0.6. 

Monte Carlo studies indicate that the B + topology consisting of two- 
prong B and three-prong D vertices has poor B + purity due to the small 
B + — > D~7i + l + v branching ratio and the large background from B° —>■ 
D(*>l + i> decays. This is corroborated by the large fraction of decays with 
B vertex charge = observed in the data for this topology. Therefore, this 
topology is rejected thereby reducing the charged sample to 634 candidates. 
These studies show that the remaining charged (neutral) sample is 97.4% 
(98.9%) pure in B hadrons with flavor contents of 66.6% B+, 22.9% B°, 
5.5% B®, and 2.4% B baryons for the charged sample, and 19.5% B + , 60.2% 
B°, 14.8% B° s , and 4.4% B baryons for the neutral sample. The sensitivity of 
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Table 1: Summary of reconstructed topologies, including the fraction of each 
topology in the combined charged and neutral sample for data and Monte 
Carlo simulation. 

B Vertex D Vertex Data MC 

# decays Fraction Fraction 

1 prong 2 prong 519 (38.0 ±1.3)% 37.6% 

Q = ±l 1 prong 4 prong 115 (8.4 ±0.8)% 8.5% 

2 prong 3 prong 149 (10.9 ± 0.8)% 9.6% 

1 prong 3 prong 341 (24.9 ± 1.2)% 26.8% 

Q = 2 prong 2 prong 175 (12.8 ± 0.9)% 13.6% 

2 prong 4 prong 68 (5.0 ± 0.6)% 3.9% 



the analysis to the individual B + and B° lifetimes can be assessed from the 
3:1 ratio of B + (B°) over B°(B + ) decays in the charged (neutral) sample. The 
fraction of misidentified leptons is 7.0% (9.7%) for charged (neutral) candi- 
dates, as determined from the simulation. However, these are predominantly 
B decays with good charge purity. 

As a check of the algorithm, the requirements on the charges of the B 
and D vertices are removed for Figs. Ua) and |l|(b). Figure |l](a) shows that, 
as expected, the charges of the lepton and D vertex are opposite for most 
reconstructed decays (provided the D vertex is charged). Furthermore, the 
charge distribution resulting from the lepton±slow transition pion vertex 
(from £)*(*)) shown in Fig. [l|(b) indicates that the track combined with the 
lepton to form a two-prong B vertex most often has charge opposite that of 
the lepton, as expected for B — > D*lv and most B — > D**lv decays. Figure [I] 
also shows the total vertex momentum distribution obtained using the tracks 
from both B and D vertices, and the D vertex multiplicity distribution (with 
the nominal charge requirements on the vertices). Overall, there is agreement 
between the data and the Monte Carlo simulation. 

The B + and B° lifetimes are extracted from the decay length distribu- 
tions of the B vertices in the charged and neutral samples (see Fig. ^) using a 
binned maximum likelihood technique. These distributions are fitted simul- 
taneously to determine the B + and B° lifetimes. For each set of parameter 
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Figure 1: Distributions of (a) the product of lepton and D- vertex charges, (b) 
sum of lepton and slow transition pion charges for data (points) and Monte 
Carlo simulation (histograms) with no charge requirement at the B and D 
vertices. Distributions of (c) total momentum of the B+D tracks and (d) 
D vertex multiplicity for data (solid circles are for total sample, open circles 
are for charged sample only) and Monte Carlo simulation (histograms are for 
total sample, shaded portions are for charged sample only). 



values, Monte Carlo decay length distributions are obtained by reweighting 
entries from generated B + and B° decays in the original Monte Carlo decay 
length distributions with W(t, r) = (^ e~*/ r J / {^— e~ t / T9en ), where r is the 
desired B + or B° lifetime, r gen = 1.55 ps, and t is the proper time of each 
decay. The fit then compares the decay length distributions from the data 
with the reweighted Monte Carlo distributions. The fit yields 



r B + 

t~b+/tbo 



1 fi~\ +{J - i6 ns 
1 56 +0 - 14 ns 

i - uo -0.14- 



with a x 2 — 78 for 76 degrees of freedom. 
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Figure 2: Decay length distributions for charged and neutral decays for data 
(points) and the Monte Carlo simulation corresponding to the best fit (his- 
tograms) . 
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Systematic uncertainties due to detector and physics modeling, as well as 
those related to the fitting procedure, are described below and summarized 
in Table |2|. A discrepancy between data and simulation in the fraction of 
tracks passing a set of quality cuts || is corrected for by removing 4% of the 
tracks from the simulation. In addition, a 0.9 m ~ 2 correction to the vertex 
reconstruction efficiency is applied, where m is the D vertex track multi- 
plicity. The uncertainty due to the track finding efficiency is conservatively 
estimated as the full difference between fits with and without these correc- 
tions. The uncertainty due to tracking resolution is similarly taken to be 



Table 2: Summary of systematic uncertainties in the B~ 
and their ratio. 



and B° lifetimes 



Systematic Error 




Ar B+ 
(ps) 


Ar B o 
(ps) 


T S 


Detector Modeling 


Tracking efficiency 
Tracking resolution 
Lepton misidentif. 




0.017 
0.020 
0.006 


0.029 
0.030 

0.007 


0.023 
0.033 
<.003 


Physics Modeling 


b fragmentation 

BR(B-*D**lv) 
BR(B-+Xlv) 

BR(B -* DDX) 

B° s fraction 

B baryon fraction 

B° s lifetime 

B baryon lifetime 

D decay multipl. 

D mom. mismatch 


0.700 ±0.011 
0.230 ±0.115 

0.15 ±0.05 
0.115 ±0.040 
0.072 ± 0.040 
1.55 ± 0.10 ps 
1.10 ± 0.08 ps 
Ref. 15 


0.035 
0.011 
0.009 
0.007 
0.008 
0.003 
<.003 
0.014 
<.003 


0.039 
0.018 
0.008 
0.007 
0.016 
0.028 
0.007 
0.009 
0.034 


0.016 
0.016 
0.011 
0.009 
0.006 
0.020 
0.005 
<.003 
0.022 


Monte Carlo and Fitting 


Fitting systematics 
MC statistics 




0.037 
0.023 


0.052 
0.024 


0.061 
0.027 


TOTAL 




0.066 


0.097 


0.088 
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the difference between fits before and after smearing and shifting the track 
impact parameters in the rz plane to account for residual VXD misalign- 
ments [^J. The smearing by a = 20 /zm/ sin# and shifts up to 20 /zm are 
required to match the core of the impact parameter distribution observed 
in the data. No correction is required to the impact parameters in the r<p 
plane. The B° lifetime is more sensitive to the above uncertainties than the 
B + lifetime because they affect the relative abundance of the various topolo- 
gies (listed in Table |l|) and the amount of wrong-charge vertices at short 
decay length is higher for two-prong than for one-prong B vertex topologies. 
The rate of lepton misidentification is varied by ±25% in the simulation. It 
was checked that the lifetimes obtained in four different azimuthal regions 
are statistically consistent. 

The 6-quark fragmentation uncertainty includes contributions from shift- 
ing the mean value of the B hadron energy |T3|| and using a different frag- 



mentation function shape [[14]]. As expected, this uncertainty affects the 
individual lifetimes but leaves the lifetime ratio relatively unaffected. Uncer- 
tainties in the B° s and B baryon production and lifetimes contribute more 
significantly to the B° lifetime, and thus affect the lifetime ratio, due to 
the larger fraction of B® and B baryons in the neutral sample. Sensitivity 
to the branching ratio for decays involving b — > c — > / transitions or for 
B — > tv t X decays is negligible. Similarly, uncertainties due to the charmed 
hadron lifetimes are negligible. 

A slight discrepancy between data and simulation is observed in the ver- 
tex total momentum distribution for the neutral sample (see Fig. |l|(c)). This 
mismatch is investigated by reweighting the Monte Carlo D vertex momen- 
tum distribution to match the data in both charged and neutral samples. 
Although the discrepancy may be attributed in part to the B decay model- 
ing, we conservatively assign an uncertainty to be the difference between fits 
with and without reweighting. 

The fitting uncertainties are estimated by varying the bin size and the 
minimum and maximum decay length cuts. Although the lifetimes obtained 
for each of these variations are statistically consistent, we conservatively as- 
sign an uncertainty equal to the root mean square value of all these results. 
This uncertainty dominates all others but is largely driven by the available 
statistics. 
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The final results are 

t b + = 1.61j#it(stat) ± 0.07(syst) ps, 
T B o = 1.56±g;it(stat) ± O.lO(syst) ps, 
t b +/t b o = 1.03ig:u(stat) ± 0.09(syst). 

These results complement those obtained with an inclusive topological tech- 
nique |n|, and are in agreement with previous measurements []2], [I7| and with 



the expectation that the B + and B° lifetimes are nearly equal. 

We thank the personnel of the SLAC accelerator department and the 
technical staffs of our collaborating institutions for their outstanding efforts. 
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